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(Dated: June 4, 2008) 

We report high precision charged-particle pair (2-D) correlation analyses in the space of A(j> 
(azimuth) and Ar/ (pseudorapidity) , for minimum bias Au + Au collisions at ^s^n ~ 200 GeV as a 
function of centrality (0-80%). The intermediate transverse momenta region chosen 0.8 < pt < 4.0 
GeV/c corresponds to an emission source size ~ 2fm obtained from HBT measurements and should 
resolve substructures at the scale of ~ 2fm. The difference and the sum of unlike-sign and like-sign 
charged pairs form Charge Dependent (CD) correlations and Charge Independent (CI) correlations 
respectively. The CD displays the initial correlation at hadronization of the opposite sign pairs 
emitted from the same space-time region as modified by further medium interactions before kinetic 
freeze-out. Our analysis of the CD correlations shows approximately jet-like structure, independent 
of centrality and is consistent with the initial correlation which is predicted by Pythia (or HIJING) 
jets. The CI correlation displays the average structure of the correlated emitting sources after kinetic 
freeze-out. For the most central bins, the Arj width of the CI correlation on the near side (A0 
around 0°) is elongated by a factor ~3 destroying the jet-like symmetry. This elongation decreases 
continually with decreasing centrality and essentially restores the jet-like symmetric structure in 
the most peripheral bins. The Pythia and HIJING event generators together with a QCD inspired 
Parton Bubble Model (PBM), which motivated this analysis, are used to compare to our data. We 
discuss the arguments for substructure, surface emission, and opacity in the central fireball region. 

PACS numbers: 25.75.Gz, 12.38.Mh 



I. INTRODUCTION 

The Search for a Quark-Gluon Plasma (QGP) 
0, 0] has been a high priority task at the Relativis- 
tic Heavy Ion Collider, RHIC % Central Au + Au 
collisions at RHIC exceed [1] the initial energy den- 
sity that is predicted by lattice Quantum Chromo- 
dynamics (QCD) to be sufficient for production of 
QGP 0. Observations of substructure have histori- 
cally played an important role in advancing scientific 
progress in nuclear and particle physics. Correla- 
tions are a powerful tool in the search for substruc- 
tures. The correlations generated by particle pairs 
(e.g. [1, S B 111) have been investigated. For cen- 
tral Au + Au collisions one might expect particle 
pair correlations to be reduced by thermalization. 
However, correlations could come from two possi- 
ble sources. One source could be detectable correla- 
tions from substructures which form on the surface 
of the fireball at kinetic freeze-out [l(| EH EH EH • 
The other is correlations from initial hard scatter- 
ings that have only been modified by interactions 
with the medium of the fireball interior but not de- 
stroyed H 1 . Various models are discussed and com- 
parisons with analyses are made in Section VI. 

In this paper we analyze the Charge Independent 
(CI) and Charge Dependent (CD) correlations in the 
two dimensional space — A</> and A77 — of charged 
pairs resulting from minimum bias Au + Au colli- 
sions at y/SNN = 200 GeV. The correlations were 
formed by particles in the intermediate p t range 0.8 
< pt < 4.0 GeV/c. The two types of charge pairs 



are the unlike-sign (US) and the like-sign (LS). The 
total correlations which are physically significant are 
CI = US + LS and CD = US - LS (Section IV and 
VI). If the background (equation 2 of Section III C) 
is subtracted from these correlations we obtain the 
corresponding signal correlations (CI signal and CD 
signal). 

We have performed the present analyses using a 
model independent method so that these intermedi- 
ate pt range correlation results can be theoretically 
analyzed in any way. However, in order to extract 
the signals we are interested in, we assume that the 
background is composed of known, expected, and in- 
strumental effects. This does introduce some model 
dependence in the signal determined by subtracting 
the background from the total correlation to the ex- 
tent of inaccuracies of the background. However, 
the robustness and characteristics of the signals ob- 
tained by this method imply that these signals are 
reasonably accurate. One should note that in com- 
paring with theoretical models (Section VI) we use 
total correlations in the model and data analysis to 
eliminate any model dependence due to the separa- 
tion of signals and background. 

In this paper we present an analysis of charged 
particle pair correlations in two dimensions — A<f> 
and A 77 — based on 13.5 million minimum bias Au + 
Au events observed in the STAR detector at ^snn 
= 200 GeV [Hj] 2 . One should note that only the 
minimum bias trigger was used. 

The paper analyzes independently the correlations 
in each of the STAR minimum bias centrality bins 



There are different types of models employed 



2 A</> = 4>\ — 4>2 and Arj = r\\ — r\2 
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from to 80% [15| and then makes comparisons with 
quantitative models and draws conclusions. The 
data analysis is very similar to that used in Ref. 

the differences will be discussed. Data cuts were 
applied to make track merging effects, HBT correla- 
tions, and Coulomb effects negligible. 

The analysis leads to a multi-term correlation 
function similar to that of equations (3 + 4 + 5) 
of Section 111 C and E of Ref. Q. This multi-term 
function fits the A-qAcj) distribution well. It includes 
terms describing correlations known to be present: 
collective flow, momentum and charge conservation, 
and instrumental effects (equation (1) Section III C). 
The sum of these terms are considered background 
defined as Bk in equation (2) (Section III A-C). 
What remains are correlations which we assume to 
represent our signals; equations (3) and (4) Section 
III E. Their robustness, characteristics and signifi- 
cance in the fits are clearly consistent with their be- 
ing signals in the unlike-sign charge pairs (US) and 
like-sign charge pairs (LS). Without the signal terms 
the fits to our data are highly rejected, but with the 
addition of the signal terms the data are well fit with 
reasonable parameterizations of the signals such as 
discussed in Ref. [!|. 

This paper is organized as follows: 

Section II describes data utilized and method of 
data analysis. Section III describes parameteriza- 
tion of the data. Section IV discusses the CI and 
CD signals and has a comparison with other data. 
Section V discusses systematic, parameter, and fit 
errors. Section VI contains a discussion and com- 
parison with models. Section VII contains Summary 
and Conclusions. The last section is an Appendix A 
which contains parameters of the fits as a function 
of centrality percentage. 



II. DATA ANALYSIS 
A. Data Utilized 

The data reported here is from STAR events taken 
at RHIC during the 2004 running period for Au + 
Au collisions at ^/snn = 200 GeV. The data were 
taken using a minimum bias trigger with the full 
STAR magnetic field (0.5 Tesla). 

The experimental arrangement is very similar to 
that described in Ref. Q. The data analysis and 
parameterization methods used previously for the 
most central (approximately 0-10%) centrality re- 
gion were employed with slight modification. The 
same parameterization was successfully used in each 
of the 9 centrality bins in this analysis. These pa- 
rameters (see Section III and Appendix A) were in- 
dependently fit in each of the nine centrality bins. 



About half the data were taken with the magnetic 
field parallel to the beam axis direction (z) and the 
other half in the reverse field direction in order to 
determine if directional biases are present. The pro- 
cedure described later in this subsection for our y 2 
analyses demonstrated there was no evidence of any 
difference in the data samples from the two field di- 
rections: the x 2 distribution of the difference of the 
two field directions in the A<f> and Ai] bins were con- 
sistent with a normal distribution expected for the 
Degrees of Freedom (DOF) . There was no evidence 
for directional biases. 

The track reconstruction for each field direction 
was done using the same reconstruction program. 
Events used in the analysis were required to have at 
least 14 primary tracks lying inside < 0.5. Tracks 
that we use are required to have at least 23 hits in the 
TPC (which for STAR eliminates split tracks), and 
have pseudo-rapidity, rj, between -1 and 1. These 
are tracks that are consistent with the criteria that 
they are produced by a Au + Au interaction. The 
surviving events totaled 7.6 million for the forward 
field and 5.9 million for the reverse field. The trans- 
verse momentum selection 0.8 < pt < 4.0 GeV/c was 
then applied. The upper p t limit of [9j was raised to 
4 GeV/c since, as discussed in [l0[ (Section 1, last 
paragraph) jet contamination is negligible if you use 
minimum bias data without a jet trigger. 

Each of the nine centrality bins was treated sepa- 
rately and fit separately. In each centrality bin the 
event records were sorted according to the z (col- 
lider axis) position of the primary vertex into ten 5 
cm wide bins along z from -25 cm to +25 cm rela- 
tive to the center of the STAR TPC. This produced 
ten files for each sign of the magnetic field in each 
centrality bin. The events for the same z bin, thus 
the same acceptance, were then merged to produce 
20 files, one for each z vertex bin, for each sign of 
the magnetic field. 

The files were analyzed in two-dimensional (2-D) 
histograms of the difference in r\ (A77), and the dif- 
ference in <f) (A</>) for all the track pairs in each event. 
Each 2-D histogram for each centrality bin had 72 
A<j) bins (5°) from -180° to 180° and 38 A77 bins 
(0.1) from -1.9 to 1.9. The sign of the difference 
variable was chosen by labeling the positive charged 
track as the first of the pair for the unlike-sign charge 
pairs, and the larger p t track as the first for the like- 
sign charge pairs. Our labeling of the order of the 
tracks in a pair allows us to range over four A</> - 
A77 quadrants, and to investigate possible asymmet- 
ric systematic errors due to geometry, magnetic field 
direction, behavior of opposite charge tracks, and p t 
dependence. Our consistently satisfactory results for 
our extensive x 2 tests of the data for these quadrants 
revealed no evidence for such effects. 
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Then we compared the Acj) - Arj data for the two 
field directions on a bin by bin basis. In the reverse 
field data, we reversed the track curvature due to the 
change in the field direction, and changed the sign of 
the z axis making the magnetic field be in the same 
direction as the positive z direction. This is done by 
reflecting along the z axis, and simultaneously re- 
flecting along the y axis. In the two dimensional Acj) 
- Arj space this transformation is equivalent to a re- 
flection in A<fi and Aij. For each pair we changed the 
sign of its Acj> and Arj in the reverse field data. We 
then calculated a x 2 based on the difference between 
the forward field and the reverse field, summing over 
the Acj) - Ar/ histograms divided by the errors added 
in quadrature. The \ 2 distribution of the difference 
of the two field directions in the Acj> and Aij bins 
was consistent with a normal distribution for ran- 
dom fluctuations of a statistical nature. Therefore 
considering the above clear justification we added 
the data for the two field directions in each central- 
ity bin. 

We compared the central z vertex bins with the 
outer z vertex bins in each individual centrality bin. 
We found no evidence of differences. Therefore we 
added the files for those 10 bins in each centrality 
bin. 



B. Analysis Method 

Separate Acj) - Arj histograms were made in each 
centrality bin for the two basic pair types unlikc- 
sign charge pairs (US) and like-sign charge pairs (LS) 
from the same-event-pairs, since their characteristics 
were different. Both histograms are needed later to 
determine the CD and the CI correlations. Similar 
histograms were made with each track paired with 
tracks from a different event (mixed-event-pairs), 
adjacent in time, from the same z vertex bin and 
thus the same acceptance. This allows use of the 
technique of dividing the histograms of the same- 
events-pairs by the histograms of the mixed-events- 
pairs which strongly suppresses instrumental effects 
such as acceptance etc., but leaves small residual ef- 
fects. 

The resulting 2-D total correlation function is de- 
fined by 

C(A0,Ar?) = S(Acj>, Ari)/M(A(j), Arj) 

where S(Acj>, Ar/) is the number of pairs at the 
corresponding values of Acj)Ar/ coming from the 
same event after having summed over all the events. 
M(A(j), Ar/) is the number of pairs at the correspond- 
ing values of Acf>Ar) coming from the mixed events 
after we have summed over all our mixed events. 



C(A</>, A77) is constructed separately for US and LS 
pairs and each is normalized to a mean of 1. 

The symmetries in the data allowed us to fold all 
four Acj) - Ar/ quadrants into the one quadrant where 
both Acj) and A77 were positive. After the cuts de- 
scribed later in this subsection, we compared the 
unfolded bins to the folded average for unlike-sign 
charge pairs and like-sign charge pairs separately. 
The folded and unfolded distributions were statisti- 
cally consistent. We searched in a number of ways to 
find asymmetries in the data via extensive \ 2 anal- 
yses and observation of fit behavior; none of any 
significance were found. By folding four quadrants 
into one we quadrupled the statistics in each bin an- 
alyzed. 

Henceforth the folded data will be used for our fits 
which increases our statistics per bin by a factor of 4. 
We used the same method of \ 2 as described in Ref. 
@ Section II C to eliminate bins which exhibited 
non-negligible track merging effects. At small Acj> 
and Ar] (i.e. small space angles) track merging ef- 
fects occur. To determine the cuts needed to reduce 
these effects to a negligible level, we varied small Acj> 
and small Arj cuts. Simultaneously the \ 2 °f an a P~ 
proximate fit to the data using equations (2) + (3) + 
(4) (see Sections III) was studied as a function of the 
bins included in the fit. With larger cuts the \ 2 be- 
haved properly until one or more of the bins affected 
by merging was included in the fit. This caused a 
huge increase in % 2 , revealing that those bin(s) were 
distorted. We confirmed by visual inspection that 
track merging clearly became important in the bins 
eliminated which caused a substantial reduction in 
track recognition efficiency. The resultant cuts also 
made the HBT and the Coulomb effects negligible. 
The required removal of bins was made by the fol- 
lowing cuts: 

For the unlike-sign charge pairs (US) 0.0 < Arj < 
0.1 and 0° < Acj) < 20°, and 0.1 < Ar, < 0.2 and 0° 
< Acf> < 10° were eliminated. 

For the like-sign charge pairs (LS) 0.0 < Ar/ < 0.2 
and 0° < Acj) < 10° were eliminated. 

The track topology differs for US and LS pairs due 
to their different curvatures in the magnetic field. 
The two tracks in the unlike-sign charge pairs curve 
in opposite directions while for the like-sign charge 
pairs the two tracks curve in the same direction. 
This makes the merging characteristics different and 
requires the cuts to be different in order to make the 
track merging effects negligible. 

The fits to the data were made over the whole 
Acj> and Ar] ranges except for the above cuts. The 
data for \Aij\ > 1.5 were cut out of the fits, since 
statistics were low and variations in efficiency are 
large. However the fits create an extrapolation for 
small Ar/ and Acj> bins where we have cut and thus 
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correct the fits for the loss of the cut-out small angle 
bins. 



III. PARAMETERIZATION OF DATA. 

We wanted to obtain a set of functions which will 
fit the data well and are interpretable to the ex- 
tent practical. We utilized parameterizations repre- 
senting known, expected physics, or attributable to 
instrumentation (or other non relevant to this anal- 
ysis) effects. Any remaining terms required to ob- 
tain good fits to the data can be considered as sig- 
nals of new physical effects. Thus signal = data - 
(known and expected) effects. For the three known 
effects (elliptic flow, residual instrumental effects, 
momentum and charge conservation) appropriate 
terms were parameterized. We then found parame- 
ters for the signal terms which are necessary in order 
to achieve a good fit to our high statistical precision 
data. 



A. Parameters Related to Elliptic Flow 

The parameters related to elliptic flow [l6| were 
represented by the usual large term 2u| cos(2A0). 
We also needed a cos(4A</>) term in order to obtain 
good x 2 fits in our intermediate p t range. 

B. Instrumental Effects 

There is a </> independent effect which we attribute 
to losses in the larger r/ tracking in the TPC and 
perhaps part of a long range correlation which is of 
no interest to this analysis and cannot be analyzed 
in detail because of the smallncss of the effect. We 
utilized mixed-event-pairs with a similar z-vertex to 
take into account these effects. Imperfections in this 
procedure leave a small depression near larger Ar/ to 
be represented in the fit by a Gaussian term labeled 
"etabump amp" and "etabump width". The width 
of this effect should be independent of the charge 
of the tracks, so we constrained it to be the same 
for like and unlike-sign charge pairs to improve the 
fit stability. We found that choosing a fixed center 
of 1.25 and a fixed width ge of 1.57 for all central- 
ity bins was adequate. However the amplitude was 
allowed to vary. Thus the functional form for this 
effect (which was treated as a background) is: 

Etabump^) = Ee~( Ar >- 1 - 25 fl 2 ^ (i) 

The residual sector effects due to the lack of read- 
out in the boundary regions (gaps between the 12 



readout sectors) clearly needed corrections for the 
2001 data utilized in Ref. [9(. However in this 
analysis the division by mixed pairs was sufficient 
to make them negligible. This was due to the fact 
that in the 2004 run, the data were corrected for the 
space charge event by event. We had attributed the 
need for residual sector corrections in the 2001 run 
to space charge effects @. 



C. Correlations associated with Momentum 
and Charge Conservation 

It is important to ensure that momentum and 
charge conservation correlation requirements are 
satisfied. For random emission of single particles 
with transverse momentum conservation globally 
imposed, a negative cos(A</>) term alone can repre- 
sent this effect since random emission of single par- 
ticles results in no correlations between them. It 
has been shown [TtJ that the cos(A</>) term alone is 
correct for conserving transverse momentum when 
there are no other correlations present. However, 
the complex correlations that occur at RHIC result 
from multiple sources which are presently not under- 
stood. It was not possible to fit our data with the 
cos(A</>) term alone. Fits were rejected by ~ 40<r 
or greater for both the unlike-sign charge pairs and 
like-sign charge pairs in the 0-50% centrality region. 
This was not surprising since random emission of 
single particles with transverse momentum conser- 
vation would not lead to the particle correlations 
observed at RHIC. Therefore we suspected that a 
more complete description of momentum and charge 
conservation was required. No one has succeeded in 
solving this complex problem in closed form even in 
the theoretical case where you observe all particles. 
It appears that a solution to this problem would re- 
quire a knowledge of all the correlations present in 
a particular analysis and this is not attainable in 
practice. Hence a reasonable approach was to try to 
solve it for the tracks we are observing in order to 
obtain a good fit. For the two variables we have, A(j) 
and A77, we used Fourier analysis and polynomial 
expansion respectively. 

Assuming that the cos(A(/>) term for random sin- 
gle particle emission was the first term in a Fourier 
expansion of odd terms, a second term cos(3A0) 
was added and found to account for almost all of 
the ~ 40<7 rejection. Based on the residual analysis 
we concluded the remaining few % required (A77) 2 
dependent terms for its removal in order to obtain 
good fits. Therefore we multiplied terms of the type 
cos(A0) and cos(3A0) by a (A7/) 2 term which when 
added reduced the remaining rejection and led to 
good fits. 
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FIG. 1: The folded - after - cuts unlike-sign charge pairs (US) correlation data vs. centrality. 




FIG. 2: The folded - after - cuts like-sign charge pairs (LS) correlation data vs. centrality. 



In addition we found that we needed a (A?7) 2 term sections A, B and C we obtain the following for our 
in our background. Such a term was added to the background Bk. 
background fit parameters since it is probably part of 
a long range correlation not relevant to this analysis. 

If we take the sum of the terms described in sub- Bk = (Known + Expected)Effects = 
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FIG. 3: "(Color online)" The signals within the prior FIG. 1-2 US (equation 3) and LS (equation 4) are the basic 
correlation data building blocks from which the two physically significant correlation signals the Charge Independent 
(CI = US + LS) signal and the Charge Dependent (CD = US - LS) signal are built (Section IV). The signals 
are obtained by subtracting the background (Bk = equation 2) from the observed total correlation. In order to 
compare different centrality bins or different experiments one must remove the dilution factor of 1/multiplicity of the 
observed correlation signals caused primarily by the quadratic increase of pair combinations. This is accomplished 
by multiplying the correlation by the multiplicity given in the last column of Table I and is defined as multiplicity 
scaling or multiplicity scaled. This figure displays a 2-D perspective plot of the multiplicity scaled unlike-sign charge 
signals (background subtracted) given in equation 3 multiplied by the multiplicity (multiplicity scaled) plotted as a 
function of A(f> and Ar/ vs. centrality. 



Boo + Bo-ii^v) 2 + B w cos Acf> + B 12 (Ai]) 2 cos A</> 
+2vj cos(2A</>) + B 30 cos(3A</>) + B 32 (Ar]) 2 cos(3A</>) 
+B 40 cos(4A</>) + Etabump(Z\77) 



D. Fitting with Bk 

We used the well known result [HI that for a 
large number of degrees of freedom (DOF), where 
the number of parameters is a small fraction of DOF 
and the statistics are high, the \ 2 distribution is nor- 
mally distributed about the DOF. The significance 
of the fit decreases by la whenever the value of \ 2 
increases by ^J2{D0F) which for our 521 DOF is 
equal to 32. 

If we fit the functional form of the background Bk 
(equation (2)) to US in each centrality bin the fits 
are rejected by ~ 50a or greater. If we fit LS to the 
background Bk in each centrality bin the fits are 
rejected by ~ 17a or greater. 



E. Signal Terms and Multiplicity Scaling 



Many signal terms in physics are Gaussian-like. 
We therefore tried fitting the signal data using two 
dimensional (2-D) Gaussian or approximate Gaus- 
sian parametric forms which successfully parame- 
terized our signals in the previous central produc- 
tion analyses [9( . The values of the signal and back- 
ground (Bk) parameters in each centrality bin were 
independently fit for that particular centrality bin. 

The physical characteristics of elliptic flow have 
been extensively investigated. They are reasonably 
understood and are essentially charge independent. 
Therefore the same large flow v 2 term was used in 
US and LS fits and its value determined as part of 
our best fit to the data set. The US in each centrality 
bin shown in Fig 1 were well fit (3<r), by adding to 
Bk an additional 2-D approximate Gaussian in Ar/ 
and A(j) given by: 
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FIG. 4: "(Color online)" a 2-D perspective plot of the multiplicity scaled like-sign charge signals (background sub- 
tracted) given in equation 4 vs. centrality plotted as a function of A<j> and A77. 



Unlike-signSignal = ^ e -(W)7^J+(^)V^-/(^) 

(3) 

Considering the enormous improvement in fit 
quality afforded by the addition of this signal term, 
we conclude that this function in equation (3) pro- 
vides a compact analytic description of the signal 
component of the unlike-sign charged pairs correla- 
tion data. The fit was improved by the addition of 
a term dependent on (A?7) 4 in the exponent (called 
"fourth" ). This was previously found to be the case 
inRef. [|. 

Our normalization is such that the correlation is 
normalized to a mean of 1. When one compares 
different ccntralitics the signal term is proportional 
to the number of correlated signal particles within 
our cuts. This is so since when two particles form 
a correlated pair they are in general not correlated 
with the remaining signal particles. Thus the pool 
of correlated particles is reduced by two every time a 
pair of correlated particles is picked. Hence the sig- 
nal term is proportional to the number of correlated 
signal pairs which is 1/2 the number of signal parti- 
cles within our cuts. On the other hand the number 
of entries to the correlation calculation grows as the 
number of particles squared. This quadratic increase 
dilutes the signal by a large factor of 1/ (particles). If 
we multiply the signal by efficiency corrected multi- 
plicity (column 3 table I) we cancel this dilution. For 



all signal comparisons as a function of centrality and 
comparisons with other experiments we utilize signal 
X multiplicity or equivalently multiplicity scaling or 
multiplicity scaled. The multiplicity used is the effi- 
ciency corrected multiplicity. The observed average 
multiplicity in the TPC and the efficiency corrected 
average multiplicity for each centrality is given in 
table I. The fits for the multiplicity scaled US signal 
data as a function of centrality for the folded - after 
- cuts data are shown in Fig. 3. 



The LS data which also could not be fit by Bk 
alone, were well described (3<r) when we added (see 
Fig. 2) a positive 2-D Gaussian and a small negative 
2-D Gaussian dip given by: 

Like - sign Signal = ^^-((^V^+^JV*^,) 

+7 4 de -((A< ! i) 2 /2<Tj d + (Ar,) 2 /2<T^) ( - 4 ) 

The significance of the small dip is ~ 20a in the 
0-30% centrality region. This dip is a physical effect 
not due to track merging. Since gluons prefer to 
emit US (pairs) and suppress emission of LS (pairs) 
in the same phase space region we attribute the dip 
in the LS as likely due to this suppression. The US 
(pairs) do not have a dip as expected. 

Therefore, we conclude that equation (4) provides 
an efficient description of the signal component of 
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FIG. 5: "(Color online)" The \ 2 for both US and LS vs. centrality using DOF as 521 and a lcr degradation of the 
fit as equal to a change of \ 2 of 32 compared to 521 results in all fits being consistent with 3cr or less. Some details 
relevant to the above procedure are given in Section III D, and Section V. 



Centrality Average Multiplicity TPC Corrected 


to 5% 


216 


292 


5 to 10% 


180 


237 


10 to 20% 


140 


176 


20 to 30% 


98 


120 


30 to 40% 


65.4 


78 


40 to 50% 


41.5 


48.2 


50 to 60% 


24.4 


27.8 


60 to 70% 


13.1 


14.6 


70 to 80% 


6.2 


6.9 



TABLE I: The average number of particles detected per 
event with 0.8 < pt < 4.0 and \rj\ < 1.0 for the 9 central- 
ity bins are given for the observed TPC particles (middle 
column), and the efficiency corrected multiplicity (last 
column). 



the LS data. The large signal is referred to as 
"lump" in the LS and is a 2-D Gaussian centered 
at the origin. It is accompanied by a small nar- 
rower 2-D Gaussian "dip" (also centered at the ori- 
gin) subtracted from it. Fig. 4 shows a 2-D A<f>Arj 
perspective plot of the multiplicity scaled LS signals 
(background subtracted) as a function of centrality. 

All the above fits are consistent with significance 
of 3cr. Note that both US and LS signal fits show a 



A?7 elongation in the central region such that the cor- 
responding angle is much larger than the A</> angle. 
These A77 elongations reduce as centrality decreases 
and virtually disappear in the peripheral region lead- 
ing to approximately jet-like symmetry. 

Appendix A contains the plots of the fitted pa- 
rameters as a function of centrality which produce 
the US and the LS fits accompanied by explanatory 
material. The signals are to some degree model de- 
pendent because of the background (Bk) subtrac- 
tion, but are robust and represent the characteris- 
tics of the underlying structures. Fig. 5 shows the 
dependence of x 2 & s a function of centrality bin for 
the fits of both the US and LS. A larger value of \ 2 
above 521 means we have exceeded the DOF. The 
difference of this x 2 divided by 32 gives the number 
of er's by which the value of x 2 departs (e.g. is in- 
creased) compared to an ideal fit (x 2 = DOF). All 
the parameters fit together, signal plus background, 
represent a model independent measure of the corre- 
lation function, which can be theoretically analyzed 
in any way chosen. 

The complex multi-dimensional x 2 surface makes 
the x 2 change non linearly with the number of er- 
ror ranges (lcr) shown on the plots in Appendix A. 
Therefore, in order to determine the significance of a 
parameter or group of parameters one must fit with- 
out them, and determine by how many a the fit has 
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CI Signal CI Signal Data 




FIG. 6: "(Color online)" a) 2-D perspective plot fit to the CI signal (equation 3 + 4) in the 0-5% centrality bin (most 
central). Note that A77 is elongated by a factor ~3 compared to HIJING or Pythia jets, which agree with Fig. 7. 
b) The CI signal data that were used in the fit. 



CI Signal CI Signal Data 




FIG. 7: "(Color online)" a) 2-D perspective plot fit to the CI signal in the 70-80% centrality bin (most peripheral). 
Note the large Arj elongation has disappeared and we have approximate jet-like symmetry (as seen in HIJING or 
Pythia jets). 

b) The CI signal data that were used in the fit. 

worsened. Then one uses the normal distribution IV. CI AND CD SIGNALS 

curve to determine the significance of the omitted 

parameter(s). A. Charge Independent (CI) Signals 



We obtain our signals by subtracting the back- 
ground (Bk) given in equation (2) from the total US 
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FIG. 8: "(Color online)" The 2-D fits for the CI signal multiplied by the multiplicity as a function of centrality. The 
CI displays the average structure of the correlated emitting sources at kinetic freeze-out (see Section IV). 



and LS correlations. If we add the US signal which 
has background subtracted to the LS signal with its 
background subtracted, we obtain the CI signal = 
US signal + LS signal. The CI signal fit and the data 
with background subtracted in the most central 0- 
5% bin that was fitted are shown as 2-D perspective 
plots in Fig. 6a-b. The CI signal displays the av- 
erage structure of the correlated emitting sources at 
kinetic freeze-out. One should note the large elon- 
gation of the corresponding A-q in the 0-5% bin was 
previously observed in Ref . @ . The CI signal fit and 
the corresponding data in the most peripheral (70- 
80%) bin that were fit are shown as 2-D perspective 
plots in Fig. 7a-b. We note the large Aij elongation 
has disappeared and we have approximate jet-like 
symmetry. 

In order to compare signals in different centrality 
bins or to compare with different analyses we must 
use multiplicity scaling. In Fig. 8 we show the mul- 
tiplicity scaled CI signals for all the nine centrality 
bins. In the 0-20% region the large Arj elongation 
mostly persists. In the 20-60% region this Arj elon- 
gation gradually decreases with the decrease in cen- 
trality. In the two most peripheral bins the elonga- 
tion is gone. 

In Fig. 9 we compare the central region CI signal 
multiplied by multiplicity for the present 2004 mini- 
mum bias Au + Au analysis with the previous 2001 
central trigger analysis. They appear quite similar 



and considering the differences in the analyses and 
the errors they agree well. 

The x 2 f° r the fits to each CI centrality bin are 
consistent with a significance of 3c or better. Fig. 10 
shows the multiplicity scaled peak CI signal vs. cen- 
trality, and the integral of the multiplicity scaled CI 
signal vs. centrality. Both have their maximum at 
the most central bin, and decrease continuously by 
a factor of approximately 3 (amplitude) and 8 (inte- 
gral) as one proceeds to their minimum values in the 
most peripheral bin. A comparison with the prior 
central trigger analysis 0] shows agreement within 
the errors in the overlap region. Also shown is a 
comparison with the results of a ^snn — 130 GeV 
analysis with a lower p t cut of 0.15 GeV/c [||. We 
attribute the observed differences, especially the dip 
at the most central CI bins, to the preponderance of 
low energy particles in that analysis. 

Fig. 11 shows the At? width as a function of cen- 
trality for the different analyses, while Fig. 12 shows 
the A(f> width for the same centralities. The present 
analysis and the central trigger prior analysis agree 
well. However in the 130 GeV analysis the A<fi width 
is decreasing with increasing centrality, while in the 
200 GeV analysis the A(f> width is essentially flat 
within the errors. However the trend shows a mod- 
est increase with centrality. We attribute the differ- 
ence in behavior to the preponderance of low energy 
particles in the 130 GeV analysis. 



13 



B. Charge Dependent (CD) Signal 

If we subtract the total LS correlation from the 
total US correlation we obtain the CD correlation. 
However, it is observed that the background (Bk) 
of the two terms are close enough in value within the 
errors to cancel each other in the subtraction. Thus 
the CD signal is essentially the same as the total CD 
correlation. The subtraction of the LS in the CD is 
considered equivalent to removal of the opposite sign 
charge pairs which are not from the same space time 
region [l||. Therefore if no further interactions 
occurred after the opposite charge pairs hadroniza- 
tion the CD would represent the initial correlation 
of those opposite sign charge pairs which are emitted 
from the same space and time region. 

The CD signal and the corresponding data in the 
most central bin (0-5%) are shown as 2-D perspec- 
tive plots in Fig. 13a-b. The large Arj elongation 
found in the CI is mostly gone in the CD most cen- 
tral bin, and we find close to jet-like symmetry. In 
the subtraction in forming the CD signals the back- 
ground terms nearly cancel out leaving small resid- 
ual background terms. Therefore to take into ac- 
count this systematic we use a procedure described 
in Section V. 

The CD signal and the corresponding data in the 
most peripheral (70-80%) bin that were fit are shown 
as 2-D perspective plots in Fig. 14a-b. Here we find 
virtually jet-like symmetry. 

Fig. 15 shows the multiplicity scaled CD signal 
vs. centrality for all 0-80% bins. All nine plots have 
the same approximately jet-like shape. 

Fig. 16 compares the multiplicity scaled CD signal 
for the present 2004 minimum bias central region 
(0-20%) with the 2001 STAR/RHIC central trigger 
analysis. They agree well. 

Fig. 17 shows the multiplicity scaled CD sig- 
nal peak amplitude vs. centrality. The CD signal 
peak amplitude shows a decreasing trend towards 
the most peripheral bins. A comparison with the 
2001 analysis (run 2) is quite good. 

Fig. 18 shows the CD Acf> width vs. centrality. 
Within the errors it is independent of centrality. The 
comparison with the 2001 central trigger analysis is 
in good agreement and the 130 GeV analysis is also 
approximately within the errors. 

Fig. 19 shows the CD Arj width vs. centrality. 
Within the errors it is independent of centrality. The 
comparison with the 2001 central trigger analysis is 
in good agreement and the 130 GeV analysis is also 
within the errors. 

All the present analysis fits for the CD quantities 
as a function of centrality are consistent with 3cr sig- 
nificance. The CD shape parameters do not change 
much with centrality (Fig. 18-19), and have close to 



jet-like behavior. 

The CD at the time of formation (hadronization) 
is a measure of the initial average distribution of the 
correlation between the US particle pairs emitted 
from the same space time region. This initial cor- 
relation distribution is expected to be jet-like and 
very similar to Pythia [H and HIJING [H| jets. If 
there were appreciable interaction with the fireball 
medium before kinetic freeze-out the initial jet-like 
correlation observed would be changed. However the 
observed CD correlation is approximately consistent 
with Pythia and HIJING jet correlations at all cen- 
tralizes (see Fig. 18-19 and 22). Therefore we can 
conclude that the emitted pairs have little further in- 
teractions after hadronization. Thus surface or near 
surface hadronization and emission from the fireball 
both occur in the central region and all other cen- 
tralizes where there is appreciable particle density. 
In the most peripheral bins the particle density is 
low enough to allow undisturbed fragmentation and 
thus no change in the CD correlation. Thus the CD 
behavior is consistent with a surface emission model 
such as Ref. [Tol j . 

V. SYSTEMATIC, PARAMETER, AND 
FITS ERRORS 

Systematic errors were minimized using cuts and 
corrections. The cuts (see Section II B) were large 
enough to make contributions from track merging, 
Coulomb, and HBT effects negligible. Systematic 
checks utilized \ 2 analyses which verified that the 
analysis results did not depend on the magnetic field 
direction, the vertex z coordinate, or the folding pro- 
cedures (see Section II ). By removing the tracks 
with |?7| >1.0 we keep the s yste matic errors of the 
track angles below about 1° [22| . 

In Section IV of Ref. @ we referred to a sim- 
ulation in Ref. [ll[ which estimated that the back- 
ground resonance contribution to the CD correlation 
cannot be more than 20% (see appendix B of Ref. 
9]). The conclusion was that the effect of these res- 
onances was primarily to increase the amplitude of 
the CD. From the resonance calculations discussed 
in appendix B of Ref. 0] the background resonance 
contribution to the US correlation was estimated to 
be 5%. Both the Parton Bubble model [l(| and HI- 
JING [2l[ only have the background resonance con- 
tribution coming from the soft beam jet fragmen- 
tation. For the present pt range model calculations 
[lol |2~H we estimated the resonance background con- 
tribution is less than 1%. 
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FIG. 9: "(Color online)" A comparison of the multiplicity scaled CI signal in the central region (0-20%) for the 
present 2004 run and the earlier 2001 run (run 2) central trigger data. 




I 

'0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Centrality 



FIG. 10: "(Color online)" The multiplicity scaled peak CI signal vs. centrality (blue open circular points). The blue 
solid circular points show the integral of the multiplicity scaled CI signal vs. centrality. The black triangular points 
show that the run 2 central trigger results agree within the errors. The red open squares points show the multiplicity 
scaled CI peak signal and the red solid square points show the multiplicity scaled CI signal integral for the 130 GeV 
analysis vs. centrality. The disagreement between the present analysis and the 130 GeV analysis, especially the dip 
at the most central, is attributed to the preponderance of low energy particles in the 130 GeV analysis. 
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FIG. 11: "(Color online)" The CI A»7 width vs. centrality for the analyses indicated. RMS(Gaussian) means that a 
Gaussian with the same RMS values as the actual fit was used to determine the width values shown. The differences 
between the present analysis and the 130 GeV analysis are attributed to the preponderance of low energy particles 
in the 130 GeV analysis. 
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FIG. 12: "(Color online)" The CI A(j> width vs. centrality. The agreement of this analysis with Run 2 (central trigger) 
is good. The disagreement between the present analysis and the 130 GeV analysis is attributed to the preponderance 
of low energy particles in the 130 GeV analysis. 
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FIG. 13: "(Color online)" a) 2-D perspective plot fit to the CD signal in the 0-5% centrality bin (most central). Note 
that the large Arj elongation found in the CI is mostly gone in the CD most central bin, and we find close to jet-like 
symmetry. 

b) The CD signal data that was used in the fit. 
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FIG. 14: "(Color online)" a) 2-D perspective plot fit to the CD signal in the 70-80% centrality bin (most peripheral). 
Note that like the CI the large A77 elongation is not present, and we find close to jet-like symmetry, 
b) The CD signal data that was used in the fit. 



17 




FIG. 15: "(Color online)" The 2-D fits for the CD signal multiplied by the multiplicity as function of centrality. 




FIG. 16: "(Color online)" Comparison of the central multiplicity scaled CD signal with the previous 2001 (run 2) 
central triggered analysis, which agree as shown in Fig. 17-19. 
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As we shall see in the next section where we dis- 
cuss model fits, the data can be well fit by assuming 
that background resonances are negligible. 

In Section VI. Systematic Errors of Ref. @ we 
concluded based on an extensive investigation for 
that paper that the increase in the CI correlation Ai] 
width would be about 2% at most. We estimate that 
in the central region of the present experiment where 
the Arj width elongation is large this approximately 
2% estimate remains the same. As the CI elongation 
becomes smaller the error estimate is adequate since 
the CI, CD, US, and LS signals approach the same 
jet-like shape which has a negligible elongation. The 
multiplicity scaling errors are estimated to be about 
2-3%. 

We have used the multi-parameter fitting proce- 
dure in the large DOF region that wc employed 
in Ref. [§]. The reader is referred to Ref. [HI 
from which we quote: "For large n(DOF), the x 2 
p.d.f. (probability density function) approaches a 
Gaussian with mean = n and variance 2 (a 2 ) = 2n." 
For n> 50-100 this result has been considered appli- 
cable, and it remains applicable and becomes more 
accurate as n increases toward infinity. It is impor- 
tant to note that the above has been proven under 
the assumption that we do not know the underlying 
physics which is certainly the case for our analyses 
and almost all heavy ion analyses at RHIC or else- 
where. In the US and the LS fits (with 521 DOF) 
an increase in \ 2 °f 32 reduces the fit significance 
by lcr (y/2DOF). The statistical significance of any 
data analysis fit in this paper can be obtained by the 
following procedure: The number of cr's of reduction 
of a fit significance = value of (x 2 - DOF)/32. The 
number of cr's refer to the normal distribution curve 
and give the probability of the fit compared to an 
ideal fit where \ 2 — DOF. For large DOF (bins - 
parameters) fluctuations occur because of the many 
bins. When one fits the parameters, they will try 
and describe some of these fluctuations. Therefore 
we need to check whether the fluctuations in the data 
sample are large enough to significantly distort the 
parameter values. This has traditionally been done 
by using the confidence level tables vs. \ 2 which al- 
lows a reasonable determination of the fluctuations 
due to binning. The approximation we have used is 
an accurate extension of the confidence level tables 
for the high (> 100) DOF region we are analyzing. 

Let us consider the method of assigning system- 
atic error ranges to the parameters. Our objective 
is to obtain systematic error ranges which are not 
likely to be exceeded if future independent data sam- 
ples taken under similar conditions are obtained by 
repetition of the experiment by STAR or others. We 
allow each parameter, one at a time, to be varied (in- 
creased and then decreased) in both directions while 



all the other parameters are free to readjust until the 
overall fit \ 2 degrades in significance by la. This 
corresponds to a change of x 2 of 32 for both the US 
and LS fits. The x 2 surface has been observed, and 
X 2 increases very non linearly with small increases of 
the parameter beyond the error ranges. Thus this is 
a conservative method of assigning systematic errors 
to the parameters determined. 

The elliptic flow has a reasonably known underly- 
ing physics which is charge independent. In the case 
of the CD elliptic flow effects have been observed to 
cancel and cause negligible uncertainties. Therefore 
we used the same large elliptic flow 2u| cos(2A0) 
term for the US and LS fits. Thus our best fit has 
this same elliptic flow term in the US and LS which 
cancels out in the CD since it is obtained by sub- 
traction of the LS from US correlation. 

The small difference of the cos(4A<ji) terms in the 
CD correlation constitute a background that must 
be determined. There are two possibilities for deter- 
mining the CD signal and background. In the first 
we defined the CD signal as the difference between 
the US and the LS signals. Thus the background is 
the difference between the US and LS backgrounds. 
This subtraction gives us amplitudes and widths de- 
termined by the % 2 surface of the US and LS fits. On 
the other hand we could directly subtract the LS cor- 
relation from the US correlation and obtain the CD 
correlation, which could also be fitted by a Gaus- 
sian plus background terms. These two procedures 
give a very similar result indicating that the small 
residual background terms are not a problem. We 
take the average of these two methods to define the 
CD signal correlation amplitudes and widths. The 
systematic errors of these measurements are deter- 
mined by exploring the x 2 surface of the direct fit 
of the CD signal correlation (method 2) plus the ex- 
tracted measurements from the difference of the LS 
signal from the US signal using the x 2 surface of the 
US and the LS fits (method 1). 

Let us now address the errors due to contamina- 
tion by including secondary particles arising from 
weak decays and the interaction of anti-protons and 
other particles in the beam pipe and material near 
the beam pipe. These secondary particle back- 
grounds have been estimated to be about 10-15% 
[23j. In this analysis we are concerned mainly with 
the angles of the secondary particles relative to the 
primaries that survive our high p t cut, not their iden- 
tity or exact momentum magnitude. Our correla- 
tions almost entirely depend on angular measure- 
ments of A(j) and Ar/. 

In the range 0.8 < pt < 4.0 GeV/c we have consid- 
ered the behavior of weakly decaying particles and 
other non-primary particles which could satisfy our 
distance of closest approach to the primary vertex 
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and p t cuts. Because high p t secondary particles are 
focused in the same direction as the primaries, only 
a fraction of these particles have sufficient change in 
angle that would cause an appreciable error in our 
correlations. We utilized the methods developed in 
Ref. [H to estimate that only 1/4 of the 10-15% sec- 
ondary particles satisfied our cuts and distance of 
closest approach. This led to our estimate that there 
is a systematic error of about 5% due to secondary 
particles which is mostly an overall scale factor. 



VI. DISCUSSION AND MODEL FITS. 

A. Discussion 

The multiplicity scaled CI signal vs. centrality is 
shown in Fig. 8. as a 2-D perspective plot. The 
CI signal displays the average structure of the corre- 
lated emitting sources at kinetic freeze-out. We an- 
ticipate kinetic freeze-out will occur at or near the 
outer surface of the fireball. The CI characteristics 
vs. centrality are shown in Fig. 10 (signal amplitude 
and integral), Fig. 11 (At; width), and Fig. 12 (A</> 
width) . 

The CI large At] elongation, signal amplitude, and 
signal integral all have their maxima in the most cen- 
tral bin (0-5%) and decrease with decreasing central- 
ity. In the most peripheral bins jet-like angular sym- 



metry is restored and HIJING [2l| can be expected 
to fit. However the central bins strongly reject a HI- 
JING fit due to the Ar/ elongation which does not 
have a jet-like angular symmetry which HIJING re- 
quires. Thus only the last few peripheral bins are 
consistent with HIJING. HIJING can not account 
for the strength of the Ar] and Acf> correlations in 
the most central collisions see Ref. [l(J. 

The CI Acj) width remains approximately constant 
with centrality (~ 30°). This implies the source we 
are viewing covers a fraction of the <fi angular range 
at all centralities in our observed pt range. This was 
calculated in Ref. [ll[ and shown to be so. In our 
Pt range the limited angular spread of particles in 
A<fi could be due to phase space focusing by flow 
as considered by Ref. @ as an explanation of the 
HBT results. A blast wave fit including all pions 
gives source sizes which are consistent with the HBT 
results for Au + Au central collisions at ^snn = 200 
GeV @. 

In the blast wave fit used in the PBM the surface 
of the blast wave is moving with maximum velocity 
(3c/4) at kinetic freeze-out when the particles are 
emitted. This phase space focusing of source size 
in the blast wave surface with the addition of bub- 
ble substructure on or near the surface of the blast 
wave is used in the PBM [lo| . This model assumes 
a bubble substructure of about a dozen similar ad- 
joining spherical bubbles in an 8 fm radius ring at 
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FIG. 18: "(Color online)" The CD Acfi width shown vs. centrality. The CD Acj> width is independent of centrality 
for 200 GeV. The present and previous analyses agree within the errors. 



mid-rapidity, perpendicular to and centered around 
the beam direction, on the surface of the fireball at 
kinetic freeze-out. The correlation functions used in 
our analysis and those employed for HBT both add 
all the bubbles in the ring on top of each other re- 
sulting in one average bubble with a source size of 
~2fm radius. If we were just dealing with flow focus- 
ing we would not observe the signal correlations we 
find in this analysis. One needs A</> regions that pro- 
duce more particles (bubble) and other A0 regions 
which produce fewer particles (background) in order 
to g enerate the observed signal correlations. In Ref. 
[10 ] a detailed QCD inspired model was developed 
which assumed the bubbles were composed of glu- 
onic hot spots which appeared in the final state of 
the fireball evolution. The PBM reasonably quanti- 
tatively fit the CI and the CD of the prior central 
trigger analysis of the approximately 0-10% central- 
ity Au + Au at snn — 200 GeV for the interme- 
diate transverse momenta 0.8< pt < 2.0 GeV/c. 



B. Model fits for the CI 

The characteristics of the present analysis dis- 
cussed above suggest that the PBM could fit the 
two most central CI (total correlation) bins and HI- 
JING could fit the most peripheral bins. It appears 
that the in between bins would require a mixture of 



the PBM, HIJING, and elliptic flow for a reasonable 
fit. Future theoretical work is expected to address 
fitting the in between bins. 

Let us first address fitting the two CI (total cor- 
relation) most central bins (0-5% and 5-10%) with 
PBM and the two most peripheral bins (60-70% and 
70-80%) with HIJING. For the two CI most central 
bins we use the Monte Carlo generated PBM events 
that were used in comparing with the central trigger 
data [9(. The entire CI correlation (signal + back- 
ground) is used for comparing the analysis results 
with the model. This eliminates any model depen- 
dence on the separation of signal from background. 
In Ref. [l(| the centrality range was 0-10%, while 
here we must separate the Monte Carlo sample into 
two pieces (0-5% and 5-10%). 

The correlation derived when using all particles 
to form pairs regardless of their charge sign is equal 
to CI/2. Since our correlations are normalized to a 
mean of 1, plotting (CI/2 - 1) provides a more easily 
understood theoretical comparison in Fig. 20. 

We need to use multiplicity scaling in order to 
compare the different centralities. This compari- 
son was straightforward for exhibiting the 2-D fits 
of the CI signal vs. centrality (see Fig. 8) but be- 
comes much harder when comparing the total scaled 
CI with models. For each centrality the CI/2 varies 
about 1.0. Since it is very hard to compare 2-D 
correlations we project into a 1-D presentation by 
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FIG. 19: "(Color online)" The CD A-q width vs. centrality. The present and previous analyses agree within errors. 



dividing 1-D A<fi projections into five Arj ranges 
0.0 < At] < 0.3, 0.3 < Art < 0.6, 0.6 < Ai] < 0.9, 
0.9 < A77 < 1.2, and 1.2 < A77 < 1.5 which cover 
our entire A77 range. 

Fig. 20 compares the increasing Arj sequence to 
the PBM fits (lines) in the two most central bins 
and to the HIJING fits (dashed lines) in the two 
most peripheral bins. The horizontal axis is Ac/) in 
degrees. 

The comparison to HIJING jets is not shown for 
the two most central bins (0-5%) and (5-10%) be- 
cause HIJING has a large disagreement with the 
data. This is due to the large observed Arj elon- 
gation (by a factor of ~3 see Section IV A) in the 
CI data which is not contained in HIJING jets. This 
contradiction of HIJING in the most central bins is 
shown in Fig. 17 of Ref. [lj3] where the centrality 
range is 0-10%. 

The away side of the scaled correlation (A0 = 
90° - 180°) is mainly V2 plus some contribution from 
cos(A0) which is close to the entire background. 
Thus the near side CI signals (A<p = 0° -90°) which 
are the signals we have extracted and analyzed in 
Section IV A should rise above the value of the away 
side. This clearly occurs in each A77 bin shown in 
Fig. 20. 

From Fig. 20 it is clear that the PBM fits for the 
0-5% and 5-10% bins agree with the CI analysis re- 
sults within an average of a few percent of the total 
correlation for all A77 ranges. HIJING fits agree with 



the analysis results for the CI 60-70% and 70-80% 
bins within an average of a few % for all A77 ranges. 
Therefore the entire 2-D analyzed CI correlations, 
which are the sum of the five Arj ranges, agree with 
the PBM fits for the 0-5% and 5-10% bins within 
an average of a few percent of the CI total correla- 
tions. The 2-D analyzed CI correlations agree with 
HIJING fits for the 60-70% and 70-80% bins within 
an average of a few percent of the CI total correla- 
tions. As we shall see in the following subsection the 
CD of the PBM agrees with the CD analysis results 
with similar precision (within a few percent of the 
correlation) . 



C. Model fits for the CD 

We again use the CD (total correlation) for com- 
parison with theoretical models. Fig. 21 shows a 
comparison of Multiplicity x CD analysis results as 
a function of A<f> (degrees) with the PBM fits (lines) 
for 0-5% and 5-10% centrality bins for the A-q range 
0.0 < A-q < 0.3. Also shown is a comparison with 
the HIJING fits (dashed lines) for the 60-70% and 
70-80% centrality bins. From the above comparison 
we observe that the PBM fits for the 0-5% and 5- 
10% centrality bins agree within a few % with the 
CD analysis results. The same is also true for the 
HIJING fits 60-70% and 70-80%. 

In Fig. 22 we plot the Multiplicity x CD as a func- 
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FIG. 20: For the Atj ranges 0.0 < Atj < 0.3, 0.3 < Atj < 0.6, 0.6 < Arr < 0.9, 0.9 < Ar; < 1.2, and 1.2 < Atj < 1.5 
we compare the multiplicity x (CI/2-1) as a function of Acj> (degrees) with the parton bubble model (PBM) fits 
(lines) for 0-5% and 5-10% centrality bins. Also a comparison with the HIJING fits (dashed lines) for 60-70% and 
70-80% centrality bins is shown. For the 5-10% the PBM fit appears to have a larger signal, but this difference is 
due to elliptic flow. The data correlations and PBM fits are shifted up by 12 units for 0.0 < Atj < 0.3, 9 units for 
0.3 < Atj < 0.6, 6 units for 0.6 < Atj < 0.9, 3 units for 0.9 < Atj < 1.2, and no shift for 1.2 < Atj < 1.5. In 
comparison with PBM in this figure we utilize (CI/2 - 1) where the CI correlation = US + LS. The comparison to 
HIJING jets is not shown for the two most central bins (0-5%) and (5-10%) because HIJING has a large contradiction 
of the data. This is due to the large observed Atj elongation (by a factor of ~3 see Section IV A) in the CI data 
which is not contained in HIJING jets. This contradiction of HIJING in the most central bins is shown in Fig. 17 of 
Ref. [Tot ] where the centrality range is 0-10%. 



tion of A<p (degrees) for 4 At] bins covering the range 
0.0 < A-q < 1.2. Each of the above four centralities 
shown was scaled so that the 5—10° Acf> bin for 
0.0 < A77 < 0.3 is normalized to 1. We did this to 
show that the CD shape is approximately indepen- 
dent of centrality. This is consistent with Fig. 18-19 
which show the Gaussian rms width in A(f> and Arj 
for the CD signal. One can see in Fig. 22 that the 
analysis points in each Arj range cluster around the 
four lines which correspond to each of the four Ar] 
ranges generated by Pythia jets [20j | . 

Pythia jets are used as the jets in HIJING and 
Pythia fragmentation is used in the PBM; thus the 
CD shape is totally given by Pythia and is approx- 



imately independent of centrality (see Fig. 18-19 
). If appreciable further interaction with the fire- 
ball interior occurred these initial CD correlations 
would be changed. Therefore we can conclude that 
the emitted pairs have little further interactions af- 
ter hadronization. Thus surface or near surface 
hadronization and emission from the fireball both 
occur in the central region and all other centrali- 
ties where there is considerable particle density. In 
the most peripheral bins the particle density is low 
enough to allow undisturbed fragmentation and thus 
no change in the CD correlation. It should be noted 
that the CD of the PBM in the 0-5% and 5-10% 
centrality bins is also reasonably consistent (within 
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a few percent of the correlation) with Pythia jets. 
Thus the CI and the CD of the PBM in the 0-5% 
and 5-10% centrality bins are in reasonable quanti- 
tative agreement with the data analysis. 

A theoretical pQCD calculation [24[ concluded 
that minijets formed with initial parton transverse 
momenta of around 3 GeV/c (also applicable down 
to 2 GeV/c) would become thermalized in a 200 
GeV/c U + U collision at RHIC and likely not escape 
from the system. However the incomplete treat- 
ment of the fireball medium in this model [24[ could 
raise some doubt about the calculation. The authors 
freely admit that they cannot use pQCD for softer 
minijet calculations. Nevertheless their pQCD cal- 
culations show that the minijets formed by initial 
hard scatterings should be thermalized using rea- 
sonable cross sections. Thus the initial correlations 
from minijets are not expected to survive to the final 
state. 

This minijet model [24[ has been modified and uti- 
lized as a qualitative picture to explain STAR cor- 



relation data by assuming that the observed corre- 
lations are surviving correlations generated by ini- 
tial hard scatterings that interact with the medium 
of the fireball interior. The correlations have only 
been modified (accounting for the change in the CI 
Ar] width see Fig. 11) and not destroyed @. This 
minijet model qualitative picture is totally different 
from the results of the pQCD treatment just dis- 
cussed which concluded the minijets are thermalized 
and do not produce the observed correlations. 

The fact that in the p t range 0.8 GeV/c< p t < 4.0 
GeV/c of this analysis the CD correlation shape is 
consistent with Pythia jets fragmentation becomes a 
serious challenge for this minijet model [![. Since the 
behavior of the CD in our p t range is consistent with 
both hadronization and kinetic freeze-out occurring 
at or near the fireball surface it is clearly implied 
that initial state hard scattering as assumed in this 
minijet model [8J] does not survive to the final state 
at freeze-out. Therefore we conclude this and likely 
any minijet model assuming survival of initial hard 
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FIG. 22: The product of the Multiplicity and the CD correlation vs. A(j> for 4 Ar; bins. Four A-q ranges with four 
centralities each are shown 0-5%, 5-10%, 60-70%, and 70-80%. Each of the centralities were scaled so that the 5 — 10° 
A<f) bin for the At) range 0.0 to 0.3 is normalized to 1. In each of the four A77 ranges we see that all 4 centrality 
bins cluster around the Pythia jets lines shown. Thus we see that the CD shape is approximately independent of 
centrality and the CD is approximately consistent with Pythia jets [2(| (dashed lines) at each centrality for all Ar/ 
ranges shown. This is true for all nine CD centrality bins (see Section IV B last paragraph), and thus is consistent 
with surface emission from the fireball at kinetic freeze-out. 



scattering to the final state correlations fails in our 
p t range. 



VII. SUMMARY AND CONCLUSIONS 

We performed an experimental investigation of 
particle-pair correlations in A<f> and Arj using the 
main Time Projection Chamber of the STAR detec- 
tor at RHIC. We investigated Au + Au collisions at 
y/SNN = 200 GeV as a function of centrality using 
13.5 million events from a minimum bias trigger to 
select the nine 0-80% centrality bins. In order to 
allow resolution of the ^2fm radii substructures we 
cut out the lower momentum particles which come 
from sources that have radii up to 6 fm in size. These 
lower momentum particles would degrade our reso- 



lution of the ~2fm sources and add particles which 
came from interacting in the interior of the fireball. 
Treatment of such interacting particles is complex. 
Thus we selected tracks having transverse momenta 
0.8 <p t < 4.0 GeV/c, and \rj\ < 1.0. We performed 
cuts to remove the small angle bins necessary to re- 
duce to a negligible level the effects of track merg- 
ing, Coulomb, and HBT. Using symmetries of the 
data in Arj and Acf> we were able to fold the four 
quadrants into one. The entire data set, unlike-sign 
charge pairs (US) and like-sign charge pairs (LS) for 
each centrality bin, was fit by a reasonably inter- 
pretable set of parameters. These consisted of 13 
parameters for the US and 15 for the LS. These pa- 
rameters are small in number compared to the total 
number of Acj> A77 bins (534 - 536) used in each cen- 
trality bin. 
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One should note that the LS signal consisted of 
two Gaussians and required only 6 signal parame- 
ters. The US signal consisted of one approximate 
Gaussian which had an additional term in the ex- 
ponent which depended on the fourth power of A77. 
Therefore four parameters were used for the US sig- 
nal. 

In order to fit the background in this precision ex- 
periment and represent known and expected physics 
effects, 9 parameters were required for both US and 
LS. The background terms are similar but not iden- 
tical in US and LS as expected. However in order 
to get the best precision and minimum distortion in 
the fits they were kept separate. The physical ar- 
guments for the necessity of these parameters were 
given in Section III. 

All the fits to the US and LS for all centrality bins 
had significance consistent with 3ct or less. The pa- 
rameterization used in this analysis was very similar 
to that successfully used in the prior central trigger 
analysis ||. The major differences were the removal 
of the readout sector gap terms and an additional 
flow term cos(4A0). Appendix A contains the pa- 
rameters and their variation with centrality. 

The Charge Independent (CI = US + LS) and 
the Charge Dependent (CD = US - LS) signals were 
formed for each of the nine centrality bins. The CI 
signal displays the average structure of the corre- 
lated emitting sources at kinetic freeze-out. Thus 
the total CI correlation is an unbiased quantitative 
measure of the particle pair correlation observed in 
the TPC. The CI signal characteristics are shown 
and discussed in Section IV A. The most central CI 
bins have the largest signals and a large Ary elon- 
gation. Both the signal and elongation decrease 
with deceasing centrality till the most peripheral 
bins where the A77 and A(f> correlation distributions 
are symmetric. The PBM fits the two most cen- 
tral bins very well within an average of a few % of 
the total correlation. HIJING fits the two most pe- 
ripheral bins with similar precision (See Section VI 
B). Future theoretical work is expected to address 
fitting the intermediate centrality region of the CI 
and to investigate the source characteristics in these 
intermediate centrality bins. 

Section V discusses systematic errors. From our 
analysis of the systematic errors we conclude that 
the conclusions drawn have not been significantly 
effected by the systematic errors. 

An overall interpretation that is consistent with 
this analysis is: 

1) The CD correlation would represent the in- 
tial correlation of opposite sign charge pairs emit- 
ted from the same space-time region if the emission 
occurred from, or very near to, the surface of the fire- 
ball at kinetic freeze-out. It was shown in Section 



VI C that the shape of the CD correlation is ap- 
proximately independent of collision centrality and 
consistent with HIJING (or Pythia) jets. Thus, we 
conclude that both hadronization and particle emis- 
sion are consistent with coming from, or very near to, 
the surface. At peripheral centralities, the interme- 
diate p t range charged particles are mainly produced 
by jets which fragment freely. As the centrality in- 
creases, a fireball forms which may be very complex, 
but the amount of fireball material that the inter- 
mediate pt particles pass through is minimal. The 
justification for this is the approximate invariance 
of the shape of the CD correlation with respect to 
centrality. Furthermore, the Pythia jet shape for the 
intermediate pt charged particles is consistent with 
the measured CD correlation for all centralities. 

2) In the most central bins, the simultaneous peak- 
ing of the 2-D CI signal amplitude, signal integral, 
and the large Arj elongation which occurs and the 
CD are all well fit by the parton bubble model. The 
PBM fits the two most central bins well within an 
average of a few percent. The model was developed 
to be consistent with the HBT observed source size. 
There is a minijet model based on pQCD calcula- 
tions for 200 GeV U + U collisions at RHIC which 
considered minijets formed with par tons of trans- 
verse momenta of around 3 GeV/c [24j . The authors 
concluded that the minijets would become thermal- 
ized and likely not escape from the system. In Sec- 
tion VI C we discussed the general characteristics of 
minijet correlations (details and limitations) origi- 
nating from hard scatterings which may be observed. 
The precision analysis of the data as a function of 
centrality presented in this paper could stimulate 
other new physics models as possible explanations 
of the observed correlations. 
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APPENDIX A: PARAMETERS 

The parameterization was similar to that em- 
ployed in Ref. 

There are two basic pairs of particles: unlikc-sign 
charge pairs (US) and like-sign charge pairs (LS). 
Both were parameterized for each centrality bin in- 
dependently. The folded after cuts US data as a 
function of centrality shown in Fig. 1 was parame- 
terized using the US parameterization. The folded 
after cuts LS data as a function of centrality shown 
in Fig. 2 was parameterized using the LS param- 
eterization. Thus there were nine independent US 
and LS sets of parameters. 

The significant difference from Ref. [9( run 2 pa- 
rameters were: 

1) The sector gap terms were not needed because 
there was event by event space charge distortion cor- 
rection for this run (run 4). In Ref. we had at- 
tributed the necessity of the sector gap corrections 
to space charge. This explains why they were not 
necessary in run 4. 

2) We found it necessary to include a cos(4A0) 
flow term in addition to the previously employed 
2v\ cos(2A0) term to obtain good fits. 

Parameters fitting requires a mixture of knowl- 
edge and experience. One must keep in mind several 
important points. 

a) One must use one's knowledge and experience 
to pick realistic parameters which will efficiently and 
meaningfully describe the data well enough so that 
the statistical significance of the final fits are good 
enough to be credible. 

b) Experienced data analysts usually require that 
their fits be consistent with at least a 3<r fit compared 
to an ideal fit in which \ 2 = degrees of freedom 
(DOF). 

1. Parameters of the background terms for US 
and LS vs. Centrality 

The first background term is an overall normal- 
ization term of the correlation to a mean of 1 shown 
in Fig. 23. 

Four terms Fig. 24-27 are needed for momentum 
and charge conservation in order to obtain accept- 
able fits. 



Two terms Fig. 28-29 represent elliptic flow ef- 
fects and are background parameters as far as this 
analysis is concerned. No reaction plane is assumed 
in this analysis. 

Fig. 30 is a background term probably due to 
remaining long range correlation. Long range corre- 
lations should be A<p independent like this term and 
important in the soft particle pt range which we do 
not explore in this analysis. 

Fig. 31 is a independent small effect which we 
attribute to losses in the larger r\ tracking in the 
TPC and perhaps part of a long range correlation 
not relevant to this analysis. We utilized mixed- 
event-pairs with a similar z-vertex to take into ac- 
count these losses. Imperfections in this procedure 
leave a small Gaussian bump near larger A77. We 
found that choosing a Gaussian with a fixed center 
of 1.25 and a fixed width of 1.57 for all centrality 
bins was adequate. 



2. Signal Parameters for the US vs. Centrality 

Fig. 32 is the US signal amplitude for the approx- 
imate Gaussian fit. 

Fig. 33 is the US signal Ar\ width for the approx- 
imate Gaussian fit. 

Fig. 34 is the US signal additional term (Ar;) 4 in 
the Gaussian exponent which makes it an approxi- 
mate Gaussian fit. 

Fig. 35 is the US signal A<p width for the approx- 
imate Gaussian fit. 



3. Signal Parameters for the LS vs. Centrality 

Fig. 36 is the LS signal amplitude for the large 
Gaussian in the fit. 

Fig. 37 is the LS signal Ar] width for the large 
Gaussian in the fit. 

Fig. 38 is the LS signal A<fi width for the large 
Gaussian in the fit. 

Fig. 39 is the LS signal amplitude for the dip 
Gaussian in the fit. 

Fig. 40 is the LS signal Arj width for the dip 
Gaussian in the fit. 

Fig. 41 is the LS signal A</> width for the dip 
Gaussian in the fit. 
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FIG. 23: "(Color online)" The normalization constant compared for US and LS. 




20 



0.004 - 



CO 



0.002 




-0.002 



-0.004 



-0.006 



_• Unlike sign pairs 
■ Like sign pairs 




i i i i i i 



i i i i i i 



i i i i i i i i i i i i i i i i 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

Centrality 



FIG. 25: "(Color online)" The cos(3A</>) term compared for US and LS. 
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FIG. 26: "(Color online)" The (A??) 2 cos(A^) term compared for US and LS. 
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FIG. 27: "(Color online)" The (Ar]) 2 cos(3A</>) term compared for US and LS. 
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FIG. 28: "(Color online)" The wa term plotted is actually 2uf and is constrained to be the same for both US and LS. 
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FIG. 29: "(Color online)" The cos(4A0) term compared for US and LS. 
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FIG. 31: "(Color online)" A small independent effect is compared for US and LS. This small term is attributed 
to a correction for the large r\ tracking effect in the TPC which remains after division by mixed pairs, and perhaps 
part of a longe range correlation not relevant to this analysis 
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FIG. 34: "(Color online)" the US signal additional term (Arj) 4 in the Gaussian exponent which makes it an approx- 
imate Gaussian fit. 
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FIG. 39: "(Color online)" The LS signal amplitude for the dip Gaussian in the fit. 
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FIG. 40: "(Color online)" The LS signal Ar/ width for the dip Gaussian in the fit. 



